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Abstract. We report on the results of a deep 1.6 Ms INTEGRAL observation of the Cassiopeia region performed from 
December 2003 to February 2004. Eleven sources were detected with the imager IBIS-ISGRI at energies above 20 keV, includ- 
ing three new hard X-ray sources. Most remarkable is the discovery of hard X-ray emission from the anomalous X-ray pulsar 
4U 0142+61, which shows emission up to ~150 keV with a very hard power-law spectrum with photon index T = 0.73 ± 0.17. 
We derived flux upper limits for energies between 0.75 MeV and 30 MeV using archival data from the Compton telescope 
COMPTEL. In order to reconcile the very hard spectrum of 4U 0142+61 measured by INTEGRAL with the COMPTEL 
upper limits, the spectrum has to bend or break between ~75 keV and ~750 keV. IE 2259+586, another anomalous X-ray 
pulsar in this region, was not detected. INTEGRAL and COMPTEL upper limits are provided. The new INTEGRAL sources 
are IGR J00370+6122 and IGR J00234+6144. IGR J00370+6122 is a new supergiant X-ray binary with an orbital period of 
15.665 ± 0.006 days, derived from RXTE All-Sky Monitor data. Archival BeppoSAX Wide-Field Camera data yielded four 
more detections. IGR J00234+6144 still requires a proper identification. Other sources for which INTEGRAL results are pre- 
sented are high-mass X-ray binaries 2S 0114+650, y Cas, RX J0146.9+6121 and 4U 2206+54, intermediate polar V709 Cas 
and 1ES 0033+595, an AGN of the BL-Lac type. For each of these sources the hard X-ray spectra are fitted with different 
models and compared with earlier published results. 

Key words. - Stars: neutron - Stars: supergiants - Stars: white dwarfs - X-rays: binaries - X-rays: individuals: 4U 0142+61, 
IE 2259+586, IGR J00370+6122, IGR J00234+6144, 2S 0114+650, y Cas, RX J0146.9+6121, 4U 2206+54, V709 Cas, 
1ES 0033+595 - Gamma rays: observations 



1. Introduction 

The INTErn ational Gamma-Ray Astrophysics Laboratory 
INTEGRAL (IWinkler et alJ 12003) is ESA's currently opera- 
tional space-based hard X-ray/soft gamma-ray telescope. Since 
its launch in October 2002 it has been observing the sky at pho- 
ton energies between 3 keV and 2 MeV. Aboard INTEGRAL 
there are two mai n hard X-ray/soft g amma-ray instruments, 
the i mager IBIS JUbertini et alJ Eo03l) and the spectrometer 
SPI ( Vedrenne et al. 2003), supplemented by two X-ray mon- 
itors, JEM-X jLund et alJl2003l) and an optical monitor, OMC 
jMas-Hesse et alJl2003h . These instruments are all wide-field 
instruments with observational coverage over a broad spec- 
tral band. In this work only INTEGR AL results from the low- 
energy detector of IBIS, called ISGRI jLebrun et alJ2003l) . will 
be presented. IBIS-ISGRI is a co ded-mask instrum ent with 
arcminute-resolution (12' FWHM: Igtos et alll2003l) imaging 
capabilities in the energy range from ~20 keV to ~300 keV. 
Its field of view (FOV) is 29° x 29° (full width zero response). 
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The most sensitive part of the FOV of a coded-mask instru- 
men t is the fully-cod ed FOV, which is 9° x 9° for IBIS-ISGRI 
(see|in 'tZandll 19921 for detailed information on coded-mask 
instruments). 

The main targets of our INTEGRAL observations are the 
young Galactic supernova remnants Cassiopeia A and Tycho. 
In this pap er these sour ces will not be addressed in detail; 
we refer to |Vink (2005) and Re naud et al.l (12004) for more 
detailed studies. This paper presents results on compact hard 
X-ray sources detected in the deep 1.6-Ms observation of the 
Cassiopeia region. In this direction of the Galactic plane the 
line of sight crosses the Perseus arm at a distance of ~3 kpc. 
One may expect that the majority of the sources are likely 
to be associated with this spiral arm. The detected sources 
are: (Anomalous X-ray Pulsar; AXP) 4U 0142+61, (high- 
mass X-ray binaries) 2S 01 14+650, y Cas, RX J0146.9+6121, 
4U 2206+54, the newly discovered IGR J00370+6122, (inter- 
mediate polar) V709 Cas, (BL-Lac) 1ES 0033+595, and the 
currently unidentified newly discovered INTEGRAL source 
IGR J00234+6141. INTEGRAL-flux upper limits are provided 
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Table 1. Summary of our INTEGRAL observations of the 
Cassiopeia region. For every observation are given the revo- 
lution number, year, begin and end time (tt, eg i„ and t en j) and the 
exposure time (t exp ). 



Rev. 


year 


date 


^begin 


date 




^end 


'exp (k s ) 


142 


2003 


Dec. 12 


15:46:50 


Dec. 


14 


21:11:53 


178.2 


143 


2003 


Dec. 15 


07:19:13 


Dec. 


17 


21:16:51 


206.8 


144 


2003 


Dec. 18 


07:15:25 


Dec. 


20 


15:00:42 


187.0 


145 


2003 


Dec. 21 


16:36:36 


Dec. 


23 


20:59:51 


173.8 


146 


2003 


Dec. 24 


18:50:50 


Dec. 


26 


20:23:19 


165.0 


147 


2003 


Dec. 27 


06:58:07 


Dec. 


29 


20:43:13 


206.8 


148 


2003 


Dec. 30 


06:51:06 


Jan. 


01 


15:10:42 


189.2 


161 


2004 


Feb. 07 


18:13:47 


Feb. 09 


18:17:17 


163.2 


162 


2004 


Feb. 10 


17:37:54 


Feb. 


12 


18:06:20 


163.2 



for some transient sources and a second AXP (IE 2259+586) 
in this field. 

For the tw o AXP sources, 0.75-30 MeV archival data from 
COMPTEL JSchonfelder et alJll993l) are reanalyzed, deliver- 
ing constraining upper limits for 4U 0142+61. In an attempt 
to reveal the nature of IGR J00370+6122 archival data from 
the RXTE All-Sky Monitor and the BeppoSAX Wide-Field 
Cameras have been studied, as well as a Target of Opportunity 
(ToO) observation with RXTE PCA. 



Table 2. List of all detected sources in the Cassiopeia re- 
gion ordered by right ascension. Most sources are High- 
Mass X-ray Binaries (HMXBs), either a Be X-ray Binary 
(BeXB) or a Supergiant X-ray Binary (SXB), with in addi- 
tion one Intermediate Polar (IP), one BL-Lac type AGN, one 
Anomalous X-ray Pulsar (AXP) and two supernova remnants 
(SNR). 



Source 


Exp. 


Flux (mCrab) 


Flux (mCrab) 


Class 




(ks) 


20-50 keV 


50-100 keV 




IGR J00234+6141 


1606 


0.72 


± 0.12 


1.4 ±0.3 


HMXB? 


Tycho 


1606 


0.43 


±0.12 


1.^0.3* 


SNR 


V709 Cas 


1606 


3.84 


±0.12 


2.8 ± 0.4 


IP 


1ES 0033+595 


1606 


0.69 


±0.13 


1.2 ±0.4* 


BL-Lac 


IGR J00370+6122 


1606 


0.78 


±0.13 


1.1 ± 0.4'' 


SXB 


IGR J00370+6122" 


489 


2.7 


±0.2 


3.6 ±0.7 




Gamma Cas 


1602 


3.8 


±0.2 


1.2 ±0.4* 


BeXB 


2S 0114+650 


1579 


7.8 


±0.2 


4.1 ±0.5 


SXB 


2S 0114+650° 


936 


12.1 


±0.2 


6.1 ±0.6 




4U 0142+61 


1152 


1.4 


±0.3 


7.0 ± 0.7 


AXP 


RX J0146.9+6121 


1133 


1.6 


±0.3 


1.5 ±0.8* 


BeXB 


4U 2206+54 


731 


2.8 


±0.5 


1.8 ± 1.2* 


BeXB 


Cassiopeia A 


1609 


3.54 


±0.13 


1.8 ±0.4 


SNR 



1 mCrab (20-50 keV)= 6.66 x lO" 6 ph cnrVkeV- 1 , 1 mCrab (50- 
100 keV)= 1.22 x 10~ 6 ph cnrVkeV - ' ; " optimized data selection; 
* possibly contaminated by artefact in map. Values to be used with 
care or used as upper limits. 



2. INTEGRAL observations and analysis 

Our observations of the Cassiopeia region were made dur- 
ing nine orbital revolutions between December 12, 2003 and 
February 12, 2004 (see TableQ]for details). INTEGRAL per- 
formed the observatio ns in a st andard rectangular (5 x 5) 
dithering mode Jjensen et alJl2003l) . The dat a provided by the 
INTE GRAL Science Data Centre (ISDC; ICourvoisier et alJ 
120031) were screened for Solar-flare events and erratic count- 
rate fluctuations due to passages through the Earth's radi- 
ation belts. After screening, the net exposure of the fully- 
coded FOV was 1.6 Ms, from 677 pointings called Science 
Windows (ScWs). The data were reduced using the Off-line 
Scientific Analysis (O SA) software package version 4.2 (see 
iGoldwurm e7alll2003[ for IBIS-ISGRI scientific data analy- 
sis). For each ScW deconvolved IBIS-ISGRI sky maps were 
created in different energy ranges. These deconvolved sky maps 
were combined into mosaic maps. 

The mosaics were examined for both known and poten- 
tial new sources. Consistency checks are necessary to eliminate 
spurious artefacts. In addition to determining the detection sig- 
nificance, the correct source shape (Gaussian in approximation) 
and the presence or absence of systematic background struc- 
tures in the maps were verified. In order to accept a source de- 
tection in this work, a significance above 3<x in at least two 
(consecutive) energy intervals was required. Source detection 
significances (^ F/dF for a background-dominated photon dis- 
tribution) were determined in the mosaics for known sources at 
their catalog positions, and for candidate new sources at the po- 
sitions in the sky maps with the highest detection significances. 



For each source and energy band the time-averaged count 
rate was determined by averaging the count rates from each 
deconvolved ScW, weighted by the variance. 

For converting subsequently count rates into flux values, 
the source spectra were normalized to the known Crab total 
spectrum (nebula and pulsar), using a 5 x 5 rectangular dither- 
ing observation of the Crab performed during revolution 102 
as reference. The source fluxes could, thus, be expressed in 
Crab units and converted in photon fluxes using the average 
Crab power-law spectrum {F(E) = aE~ r ) with photon in- 
dex T = 2.108 and no rmalization a = 9.59 ph ern es ' at 
1 keV JWiliingaleetalJl200ll) . This relation is valid through- 
out the IBIS-ISGRI waveband up to ~180 keV, as was verified 
bv iKuioer etail d2006l) . This method for deriving IBIS-ISGRI 
source spectra enables us to determine the spectra without de- 
tailed knowledge of the instrument response, which contained 
significant uncertainties at the time of our analyses. 

Spectra l-model fitting was performed with XSPEC ver- 
sion 11.3.1 jArnaudl 19961) . For power-law fitting the PEGPWR 
model was used, which is a normal power law, but the program 
returns the power-law index and the unabsorbed power-law flux 
(^poweriaw) m a cnosen energy band instead of a normalization 
at a given energy to describe the spectrum. This model has been 
chosen in order to minimize the covariance between the param- 
eters, resulting in the best error estimates. 

All errors quoted in this paper are lcr errors, unless stated 
otherwise. 
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Fig. 1. Cassiopeia-region significance map in Galactic coordinates (I, b) as seen with IBIS-ISGRI in the 20-50 keV energy range. 
All excesses with positive significances are drawn. Indicated with white circles and annotation are the sources significantly 
detected during these observations. Indicated with yellow circles and annotation are transient sources which have been weakly 
detected in this work and/or detected with INTEGRAL at other epochs. AXP IE 2259+586, not yet detected, is a good candidate 
for future study (see section l4~2l . The weakest (marginally) detected source is Tycho with 3.5<x, while the brightest source in this 
field, 2S 01 14+650, is detected with 43.1cr 



3. Sky maps and source list 

Figures^andlJlshow sky maps of a large part of the Cassiopeia 
region in the energy bands 20-50 keV and 50-100 keV, re- 
spectively. For presentation purposes the sky regions at higher 
Galactic latitudes and below Galactic longitude / ~ 108° 
were cropped off. As a result 4U 2206+54, located at (/, b) = 
(100?6, -1°1), is not shown in Figures^and|2] The significantly 
detected sources that passed our consistency checks as men- 
tioned in Sect. |3 are indicated by white circles and annotation. 
The yellow circles and annotation indicate weakly detected 
sources, transient sources, that were detected by INTEGRAL 
at other epochs and the second AXP in the field, IE 2259+586, 
for which there is so far no INTEGRAL detection. The weak 
(one energy interval) detections and upper limits of the lat- 
ter sources are given and discussed in Sects. 14.21 and [6] All 
(weakly) detected sources are listed in Tabled 

4. Anomalous X-ray pulsars 

Anomalous X-ray pulsars, originally dis cussed as a 
subset of low-mass X- ray binaries (e.g. iHellierl 1 19941 
iMereghetti & Stellalll995l) . are members of a special subset 
of isolated neutron stars fsee lWoods & Thompsorj2005l for a 
recent review). So far, there are seven established AXPs and 
one AXP candidate. Their spin periods (P) are found in a nar- 
row range between 5.5 s and 11.8s and their period derivatives 
(P) are of the order of ~10~ n s s _1 . Their characteristic ages 
(t ~ P/2P) indicate that they are young objects, though the 



characteristic age should be interpreted with care, especially 
for young neutron stars. However, the facts that all the Galactic 
AXPs (one is located in the SMC) are found within ~1° of 
the Galactic plane and that three of them are associated with 
supernova remnants (Gaensler et al. 2001) support their young 
nature. These pulsars are called anomalous, because their 
X-ray luminosities below 10 keV are orders of magnitude 
higher than the rotational energy losses. The magnetic field 
strengths of AXPs, B oc yP~P, are of the order 10 14 G, much 
stronger than estimated for normal (radio) pulsars (B ~ 10 12 
G). 

In order to explai n the energy budget of the A XPs, the 
magnetar model ([Thompso n & Duncan! 1 19931 1 19961) is cur- 
rently the most accepted model. This model was originally 
aimed at explaining phenomena observed in Soft Gamma- 
ray Repeaters (SGRs), but due to recent observations of sim- 
ilar behaviour of SG Rs and AXPs , like bursting activity 
dGavriil et alJl2002t iKaspi etalj Eo03i. they are now believed 
to originate from the same class of objects. This model also 
explains the gl itching p henomen a detected in the spin down 
of two AXPs (lKaspi&Gavrii]ll2003l balFOsso etaii r2003: 
IWoods et al.l2004l) . The magnetar model postulates that AXPs 
and SGRs are powered by magnetic-field decay. 

In the pre-INTEGRAL era AXPs were known as soft X- 
ray sources visible up to 10 keV. All AXP spectra, below 10 
keV, could be described by a black-body component and a 
soft power law. However, this view has changed completely 
since INTEGRAL discovered hard X-rays above 20 keV from 
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Fig. 2. Cassiopeia region presented in the same format as in Fig.^i n the 50-100 keV energy range. All detected sources are 
indicated with white circles and annotation. Several sources are no longer apparent in this energy range. IGR J00370+6122 is not 
visible in this total-exposure map, but is detected during two single revolutions (see Sect. l5.4t . Some residual artefacts/systematics 
can be seen in the center of this map. 



the dire ction of AXPs IE 1841-045 and 1RXS J 170849- 
400910 (Mol kov et alJl2Q04 iRevnivtsev et al.lEo04al) . In ad- 
dition, ^^^to^^[^04^repoitedtiielNTEGRAL de- 
tection of 4U 0142 +61, in a pre liminary analysis of the data 
used in this work. Kuip er et alJ {2004) provided definite ev- 
idence that the hard X-rays from the direction of IE 1841- 
045 indeed originated from the AXP rather than the associated 
supernova remnant G27. 4+0.0 (Kes 73). They showed that a 
significant fraction of the hard X-ray emission is pulsed, us- 
ing RXTE PCA and HEXTE covering the same energy range 
as INTEGR AL. Also for AXPs 1RXS J 170849-4009 10 and 
4U 0142+61 iKuiper et alJ §006) subsequently found pulsed 
emission in the hard X-ray regime above 20 keV using archival 
RXTE data. 



4.1. 4U 0142+61 



4U 0142+ 61 was d etected by t he Uhuru satellite in the early 
seventies (Gia cconi et alJll972t iForman et al1ll978h . The na- 
ture of 4U 0142+61 was unclear for more than two decades, 
caused by source confusion with the nearby located high-mass 
X-ray binary RX J0146.9+6121 (24' from the AXP, see sec- 
tionES). With th e discovery of RX J 0146.9+6121 in ROSAT 
0.1-2.4 keV data jMotch et alJll99ll) . lMereghetfi et alJ Jl993h 
realized that t he 25-minute peri odicity, up to then attributed to 
4U 0142+61 JWhite et al] 19871) . was unlikely to originate from 
4U 0142+61, b ecause of the characteristics of the modulations. 
iHellieJ i 1994 showed with ROSAT PSPC observations that llsrael et al] Jl994 . They found that the rotational energy loss 
this periodicity indeed originated from RX J0146. 9+6121. The of 4U 0142+61 is not sufficient to explain its X-ray luminosity 
spin period of 8.7 seconds of 4U 0142+61 was discovered by below 10 keV. 





W'\ v >r 1 


T 

... .i i . . . 


1 1 



1 10 100 1000 10' 



Energy [keV] 

Fig. 3. The broad-band X-ray to gamma-ray spectrum of 
4U 0142+61 in an E 2 F r epresentation. Shown are the 1-8 keV 
Chandra fit JPatel etalJliool . the 20-150 keV IBIS-ISGRI 
flux values (lcr error bars) with 2 fits (see Table [3} and the 
0.75-30 MeV COMPTEL 2cr upper limits. The dotted lines 
are extrapolations of the models fitted to the IBIS-ISGRI flux 
values. The extrapolation of the single power law shows that 
the spectrum has to bend/break in order not to be in conflict 
with the COMPTEL upper limits. 
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Fig. 4. Fit residuals for the 4U 0142+61 spectrum: a) for a sin- 
gle power-law fit; b) for a power-law with a high-energy cutoff 
fit. 



4U 0142+61 was the first AXP for which optical emission 
was discovered (Hulleman et al. 2000), which later turned out 
to be pulsed jKern & Martinl2002l) . lHulleman et alJ (|2004) also 
disc overed near-IR emission from this source, further studied 
bv lMoriietaT]j2005l) in simultaneous X-ray, near-IR and opti- 
cal observations. They discovered in Subaru observat ions, that 
also the near-IR emission is pulsed jMorii et alJl2006l private 
communication). So far, no radio emission has been detected 
jGaensleret alJ200lh . 

The soft X-r ay spectrum observed with Chandra by 
IPatel et alJ d2003h shows a typical AXP spectrum. The spec- 
trum can be fitted with a thermal black-body model with a tem- 
perature kT = 0.470 + 0.008 keV and a very soft power-law 
component with a photon index F = 3.40 + 0.06. 

In this work no INTEGRAL-timing analysis has been per- 
formed for4U 0142+61. This will be addressed in future work, 
exploiting also additional INTEGRAL observations. 

4.1 .1 . 4U 0142+61 : INTEGRAL results 

4U 0142+61 was significantly detected by INTEGRAL as a 
hard X-ray source (see Figs . [Hand El for a preliminary report 
see lden Hartog et al.l2004al) . even though this source is located 
17°7 from Cassiopeia A and 9? 5 from Tycho, therefore being 
always outside the fully-coded FOV of IBIS-ISGRI. In the sky- 
image mosaic with a total exposure of 1.15 Ms, 4U 0142+61 is 
detected with 5.0<x, 9.5<x and 4.6<x significances in the energy 
bands 20-50 keV, 50-100 keV and 100-150 keV, respectively. 

The IBIS-ISGRI spectrum for energies above 20 keV, 
shown in Fig. [3] can be described by a power-law model with a 
photon index F = 0.73 + 0.17. This fit has & x 2 r = 2.01 (dof = 7) 
and thus a null-hypothesis probability of 5% (see Tableland 
Fig-Et for the fit residuals). Striking is the remarkable harden- 
ing of the spectrum above 10 keV when compared to the very- 
soft Chandra spectrum, also shown in Fig. [5] The transition 
occurs between the Chandra and IBIS-IS GRI energy w indows 
where it can be studied with RXTE ( Kui per et al.1 .2006). 



4.1.2. 41) 0142+61 : COMPTEL archival-data analysis 

The Compton telescope COMPTEL dSchonfelder et alJll993h 
aboard the Compton Gamma-Ray Observatory (CGRO, 1991— 
2000) is sensitive to gamma-ray photons between 0.75 MeV 
and 30 MeV, thereby covering the harder gamma-ray band ad- 
jacent to the INTEGRAL one. The very-hard spectrum mea- 
sured with IBIS-ISGRI for 4U 0142+61 prompted us to rean- 
alyze the COMPTEL data archive to search for a signal from 
4U 0142+61. COMPTEL has a wide circular field of view cov- 
ering ~ 1 steradian. Therefore, during the long mission lifetime 
of CGRO, most of the sky, particularly the Galactic plane, has 
been viewed with long exposures. Sky maps and source pa- 
rameters can be obt ained with the maximum-likelihood method 
JdeBoeretaT]ll992l) . which is implemented in the standard 
COMPTEL data-analysis package. For the analysis and data se- 
lections the approach described by Zhang et al. ( 2004) was fol- 
lowed. The effective on-axis exposure time on this AXP, accu- 
mulated over the total mission, amounts to ~4.2 Ms. Standard 
energy intervals were selected for the analysis. In the lower- 
energy intervals, 0.75-3 MeV and 3-10 MeV, weak, low sig- 
nificance (^1.5cr) excesses were found. In the 10-30 MeV 
energy band no positive-flux measurement was derived. The 
2cr upper limits for these energy intervals are 3.1 x 10 -8 
phcirrV'keV- 1 , 2.9 x 10~ 9 ph cnrV'keV- 1 and2.0xlO- I() 
ph cm _2 s _1 keV _1 , respectively. These are shown in Fig. 3, 
and can be directly compared with the hard X-ray spectrum. 
If the power-law spectrum fitted to the IBIS-ISGRI data ex- 
tends to the MeV energy range, it would result in a photon 
flux of -6.7 x 10~ 7 ph cnrV^eV- 1 (0.75-3 MeV), i.e. a 
bright COMPTEL source, inconsistent with the COMPTEL up- 
per limits. It shows that the spectrum has to bend or break be- 
low ~0.75 MeV. To search for possible indications of this break 
in the IBIS-ISGRI data, a power-law model with a high-energy 
cutoff was fitted to the data. This resulted in a fit improvement 
of 2.34<t with a A X 2 of 7.9 for 2 extra free parameters (see 
Table and Fig. [3}. This fit gives an extremely hard power 
law in the lower energy part with a photon index F = - 0. 1 1 + 
0.43 and a high-energy cutoff at 73 ± 15 keV with an e-folding 
energy of 37 + 18 keV. With this model no significant flux is 
expected in the COMPTEL energy range. Figure|4j> shows the 
residuals from the fits. We merely consider this 2.3<x fit im- 
provement as a hint for a break. With recently performed ad- 
ditional INTEGRAL observations we expect to obtain tighter 
constraints. 

4.2. 1E 2259+586 

IE 2259+586 has been very important for AXP research. Not 
only was it the fir st of the currently known A XPs recognized as 
a peculiar source {fahlman & Gregory! 1 98 ll) . but it also played 
a key role in unifying AXPs and SGRs as magnetar candidates 
when this source s howed more than 80 SGR-like bursts in 2002 
dKaspi et all2003h . The soft X-ray spectrum, which can be de- 
scribed with a black body with a temperature kT = 0.41 keV 
and a very soft power- law with photon ind ex T = 3.6, is a 
typical AXP spectrum (see Patel et al. 2001, for recent X-ray 
observations). Interestingly, Woo ds et alJ (|2004) observed sig- 
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Table 3. Spectral fit results for AXP 4U 0142+61. The data 
can be fitted with a single power-law model. Using a power- 
law with a high-energy cutoff model, the fit improves 2.34<x, 
indicating a hint for the expected break in the spectrum. 



Power law (20-150 keV) 
T 0.73 ±0.17 

Fpoweriaw (9.7 ± 0.9) x 10-" erg cm-V 1 

Xl (dof) 2.01 (7) 

Power law with High-Energy cutoff (20-150 keV) 
T -0.11 ±0.43 

^Poweriaw (18 ± 5) x 10"" erg cm~ 2 s _1 
Cutoff 73 ± 15 keV 
Efoid 37±18keV 
Xl (dof) 1.23 (5) 

gunatg 9.5 + °; x 1Q-" erg cm- 2 s-' 




Energy [keV] 



Fig. 5. Soft X-ray spectral fits and hard X-ray to gamma-ray 
upper limits of IE 2259+586 in a E 2 F representation. Shown 
are fits measured by Chandra (0.5-7 k eVIPatel et al : 2001 ) and 
by XMM-Newton (2-10 keV; Obs. l. lWoods et aljl20Q4t) : the 
IBIS-ISGRI 20-150 keV and COMPTEL 0.75-30 MeV 2cr up- 
per limits. The solid staircase line is the IBIS-ISGRI 3cr sen- 
sitivity for 5-Ms exposure time. A hard high-energy spectral 
component as found for 4U 0142+61 is not in conflict with the 
INTEGRAL upper limits. 

nificant spectral changes (for energies below 10 keV) before 
and after the 2002 outburst. 

Contrary to 4U 0142+61, IE 2259+586 is associated 
with a young supernov a remnant G109. 1-1.0 (C TB 109), 
which was discove r ed b y Greaorv & Fahlman ( 1980) (see also 
lRho& Petrel 1 1997t IPatel et al]l200lt ISasaki et al.ll2004h . The 
supernova remnant does not influence the INTEGRAL analy- 
sis as it is only detected up to ~4 keV JSasaki et al.l2004h . 

4.2.1 . 1 E 2259+586: INTEGRAL results 

IE 2259+586 is located only 2°.9 from Cassiopeia A (see 
Fig. Therefore it was optimally exposed for 1.6 Ms in 
the fully-coded FOV. Even with the current high exposure, 
the sensitivity of INTEGRAL is not sufficient to detect this 




10 20 30 40 50 60 
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Fig. 6. IBIS-ISGRI (20-50 keV) light curve of 2S 0114+650 
with 0.1 day binning. It shows periods of flaring as well 
as quiescent phases. The epochs of maximum flux, using an 
ephemeris based on RXTE-ASM data, are indicated with ar- 
rows. The four time intervals that are omitted for the spectral 
analysis are indicated as thick horizontal bars. 



AXP above 20 keV. Only upper limits could be extracted from 
the IBIS-ISGRI data. In the energy bands 20-50 keV, 50-100 
keV and 100-150 keV the 2cr upper limits are 3.2 x 10~ 6 
ph ctrT VkeV- 1 , 1 .4 x 10~ 6 ph cirT V'keV- 1 and 0.9 x 10~ 6 
ph cm _2 s _1 keV _1 , respectively. Given the reported variability 
below 10 keV, the data were also analyzed in separate revolu- 
tions, but IE 2259+586 was not detected. 

4.2.2. 1 E 2259+586: COMPTEL archival-data analysis 

Also for IE 2259+586, analysis of COMPTEL ar chival data 
is performed as described for 4U 0142+61 in Sect. 14. 1 .21 The 
effective on-axis exposure time on IE 2259+586, accumulated 
over the total mission, amounts to ~4.9 Ms. As for4U 0142+61 
COMPTEL was not sensitive enough to significantly detect this 
AXP. The 2a- upper limits for the 0.75-3 MeV, 3-10 MeV and 
10-30 MeV bands are 2.3 x 10~ 8 ph cirrV'keV- 1 , 1.7 x 10~ 9 
ph cm _2 s _I keV _1 and 2.0 x 10 ph cm" 2 s _1 keV _1 , respec- 
tively. The 2cr upper limits are shown in Fig.|5] together with 
the IBIS-ISGRI upper limits and spectral fit s to two different 
spectra measured below 10 keV by Chandra JPatel etalJ200ll) 
and by XMM-Newton JWoods et alJl2004 . A spectral shape 
including a hard spectral tail like measured for 4U 0142+61 
(Fig. is allowed by all upper limits above 10 keV. 

Recently Kuip er et aD (|2006) detected pulsed emission up 
to -24 keV with the RXTE-PCA, indicating the likely onset 
of a hard tail, which is still consistent with the upper limits 
derived in this work. This detection is promising for a future 
INTEGRAL detection. The total exposure on IE 2259+586 
will be more than 5 Ms by the end of 2006. 
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5. Binaries 

Since launch, INTEGRAL has discovered 107 new hard X- 
ray sources 1 , among which a la rge number o f binar i es (se e 
e.g. iRevniytsev et alJ l2004albl) : iMolkov et all d2004 1200.4 : 
Bird et al J 11200 4. 2006) and references therein). One of the im- 
portant discoveries of INTEGRAL is the discovery of a pre- 
viously missed class of objects, namely highly absorbed High- 
Mass X-ray Binaries (HMXBs) (N u > 10 23 cm _1 ; see lKuulkersl 
2005, for a recent review). 

HMXBs can be divided in two classes dCorbetlll986T) : the 
Be X-ray binaries (BeXBs) and the supergiant X-ray (SXBs) 
binaries, which have an OB supergiant or a bright giant as 
donor star. Currently there are about 150 HMXBs known, 
out of which ~65 are located in our own Galaxy. Of these 
150 systems, about 20 are of the supergiant class. From the 
new INTEGRAL sources, for which follow-up observations 
are performed, already 10 are identified as supergiant X-ray 
binaries with a high intrinsic absorption. Including the new 
INTEGRAL SXBs, the s ample of SXBs has grown to a bout 
30. See e.g. iNegueruelal J2005I) : iNegueruela et alJ J2006h and 
iLutovinov et alJ (I2006 ) ) for the latest developments and impli- 
cations of the newly found HMXBs by INTEGRAL. 

In this section the binaries that were detected in the 
Cassiopeia field (see Tabled are introduced and the results are 
presented. It concerns four HMXBs, one Intermediate-Polar 
(IP) and the new source IGR J00234+6141 which is most likely 
an HMXB. 



5.1. 2S 0114+650 

2S 01 14+650 is a supergiant X-ray binar y. The supergiant com- 
panion LSI +65°010 is of the Blla type jReig et all 19961) . This 
system exhibits X-r a y modulations with a period of 2.7 hours 
jFinlev et alJ EM ICorbet et alJ EM iBonning fc Falanga 



2005) which are interpreted as originating from the neutron 
star's spin - th e slowest known to d ate. The orbital period 
of 11.59 days ( Cramp ton et al J 1 1 9851) is quite normal for a 
SXB. Recentlv. lFan-ell et alJ J2004) found evidence for a super- 
orbital period of 30.7 days, which they interpret as an indica- 
tion for the presence of a precessing warped accretion disk. 

2S 01 14+650 has be en observed in the soft X-ray band sev - 
eral times (<10keV; e.g. lHall et alJ2000tlApparao et al.ll99ll) . 
However, before INTEGRAL little informatio n was avail able 
for the high-energy part of the spectrum. iKoenigsberger et all 
ill983h published an OSO-8 spectrum taken in 1976, which 
could be fitted with a hard power law with photon index F — 1 .2 
and a high-energy cutoff E cut ~14 keV. The source was de- 
tected up to ~30 keV. 

5.1 .1 . 2S 01 14+650: INTEGRAL results 

INTEGRAL observed 2S 0114+650 for 1.58 Ms. It was de- 
tected at significance levels of 43.1cr and 8.6cr in the 20- 
50 keV and 50-100 keV bands, respectively. The light curve, 
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Fig. 7. Time averaged IBIS-ISGRI 20-150 keV spectrum of 
2S 0114+650 in an E 2 F representation. This spectrum shows 
a high-energy tail which was not observed before INTEGRAL. 
The best fit to the data is a bremsstrahlung with a power-law 
model (see Table 0}. The best-fit model and the separate com- 
ponents are plotted. 
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1 see http://isdc.unige.ch/~rodrigue/html/igrsources.html for the 
complete IGR-sources list 



Fig. 8. Residuals for IBIS-ISGRI 20-150 keV spectral fits 
of 2S 0114+650. Panel a) single power-law model; b) sin- 
gle bremsstrahlung model; c) bremsstrahlung plus power-law 
model. 



presented in Fig. [6] shows a variable source with periods of 
flaring and quiescence, both lasting of the order of less than 
a day to days. This source is known to show flaring activ- 
ity adding up to a single broad pulse w hen folding with the 
orbital ephemeris over a long time span (IKoenigsberger et alJ 
ll983HApparao et alJl99lUCorbet et al.ll999HHall et al.l200fi) . 
Using an ephemeris (maximum flux expected at MJD = 
51825^2 ± 0?3 + N x 11*5995(53)) created with the available 
online RXTE-ASM data (MJD 50087-53100), we find that the 
phases of the flares are distributed consistently over the phase 
interval of the broad time-averaged pulse. 

For the spectral analysis, four time intervals in which the 
source was in a weak/quiet state are omitted. The omitted inter- 
vals in MJD are: 52988-52990, 52992-52996, 52997-53000 
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and 53004-53006 (these intervals are indicated in Fig.[6ji. The 
remaining data consist of 374 ScWs adding up to an exposure 
of 935.8 ks. For this selection of observation times the detec- 
tion significances increased to 52.3cr and 10. Oct in the 20-50 
keV and 50-100 keV bands, respectively. 

The time-averaged high-energy spectrum (20-150 keV) of 
2S 0114+650 is shown in Fig. A fit with a single power- 
law model (F = 3.01 + 0.06) gives an acceptable fit up to 
150 keV (see Tableland Fig.|8jj). A more physical thermal- 
bremsstrahlung model could not fit the data satisfactorily up 
to 150 keV. Above 70 keV it underestimates the high-energy 
flux and the fit is therefore not acceptable (see Table |4] and 
Fig.[3J>). Adding a power-law component to the bremsstrahlung 
model, a good fit is achieved. However, statistically it can not 
be favoured over, nor distinguished from the single power-law 
model (see Tableland Fig.G^c). To explain the hardest part 
of the spectrum c omptonization models are invoked. The mod- 
els by iTitarchukl i ll 9941 compTT) and Po utanen & Svenssonl 
(1996, compPS) are able to model the high-energy part of 
the spectrum better than the bremsstrahlung model, but the 
parameters are poorly constrained. Therefore the high-energy 
tail which has been observed by INTEGRAL for the first time 
can be described by either comptonization in a thermal or non- 
thermal plasma. It is not clear whether it originates in a hot 
corona or in a putative accretion disk. The accretion-disk sce- 
nario, however, is less likely for SXB systems. 

Re cently iFihppova et alJ J2005h and IBonning & Fa langa 
( 120051) have independently obser ved 2S 0114+65 with 
INTEGRAL. Filippov a et al.l J2005I) caught the source in a 
low state and measured a simple power law (r = 2.3 ± 
0.4) which is marg i nally consistent (2cr) with our findings. 
IBonning & Falanga! J2005t) observed 2S 0114+6 for 180 ks 
in December 20 04. They find similar results compared to 
lHall et alJ feOOOh combining JEM-X and IBIS-ISGRI data. 
Their best-fit model is a power law with a high-energy cutoff: 
photon index F = 1.6 + 0.5 and high-energy cutoff E cut = 9 ± 1 1 
keV with an e-folding energy E foU = 22. l^ 1 keV (90% er- 
rors). Fitting our data with these parameters, leaving the nor- 
malization free, gives an unacceptable fit. This model also fails 
to fit the hard tail and therefore underestimates the flux above 
70 keV. Already in their spectrum it can be seen that the model 
is below their data points above 70 keV, but just because of the 
larger statistical errors their mod el was sti l l acc eptable to de- 
scribe their spectrum. IBonni ng & Falanga (2005) note the in- 
dication for a possible cyclotron absorption feature at 44 keV. 
Our data do not confirm this indication, despite the much longer 
exposure time. 

Even more recently, Mas ettietal1 l2005) published results 
from a BeppoSAX observation of 63 ks with the narrow field 
instruments, sensitive between 0. 1 keV and 300 keV. The spec- 
trum of 2S 0114+65 is determined in the 1.5-100 keV band. 
Their findings agree with the findings of IBonning & Fal anga 
J2005I) . except that they also do not find evidence for a cy- 
clotron absorption line. Compared to our spectrum their fit re- 
sults also fail to detect the hard tail. 



Table 4. Spectral fit results for 2S 01 14+650. 
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5.2. y Cas 

y Cas i s a bright B0.5 IVe star and represents the prototypical 
Be sta r JSecchi 18671). It was di scovered in X-rays bv ljerniganl 
Jl976l) and lMasonetall ( ll976l) . The X-ray characteristics of 
y Cas made it a peculiar star because its X-ray luminosity is 
several times higher than measured for other Be stars, but ~2 
times weaker than the weakest BeXB (Murak ami et alJ fT986). 
Through radia l velocity analysis of sp ectral lines, it was found 
to be a binary dHarmanec et al. 2000). 

Determining whether the nature of the compact companion 
is a neutron star or a white dwarf was a c hallenge for a long 
time. Rappap ort & van den Heuvel {l982) argued that the X- 
ray luminosity of ~10 33 erg s~' combined with a wide orbit, 
and the hardness of its spectrum, support a neutron-star com- 
panion. However, most authors agree that a white dwarf com- 
panion is more likely: its spectrum can be fitted with a ther- 
mal optically thin plasma with an extremely high temperature 
of 10-12 keV; a strong Fe emission li ne at 6.8 keV with an 
equivalent width o f 280 eV is present (Murak ami et alj Fl986: 
lOwens et al.l [T999l: and also rapid fluctuations on timescales 
of sec onds hinting towa rds a white-dwarf companion are de- 
tected JKuho et al Jl 998USmirh et al Jl 998h. 

lOwens et all £l999) observed y Cas with BeppoSAX. The 
spectrum of y Cas could be extracted up to ~30 keV with the 
HPGSPC and -40 keV with the PDS. The best fit to their X-ray 
spectrum from 0.2 -^-0 keV is achieved with a MEKAL model 
jMewe et all 19 95) with an additional Carbon line (at 0.3 keV), 
or a two component MEKAL model. The temperature of the 
hot plasma in both models is 12.4 + 0.6 keV. 

5.2.1 . y Cas: INTEGRAL results 

INTEGRAL detected y Cas with high significance up to 50 keV 
(25.1cr). Between 50 keV and 100 keV a 3<x flux point is mea- 
sured (see Table |2j, however, this flux point is likely contami- 
nated by surrounding artefacts. For this reason it was discarded 
in the spectrum, shown in Fig.|9] 
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Fig. 9. IBIS-ISGRI spectrum of y Cas in an E 2 F representation. 
Overplotted is the best-fit model, which is a power-law with a 
high-energy cutoff (see Table [3}. The bottom panel shows the 
fit residuals. 

Table 5. Spectral fit results for y Cassiopeia. 
Power law with High-Energy cutoff (20-50 keV) 
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The INTEGRAL spectrum of y Cas can best be described 
by a power-law model with a high-energy cutoff. However, the 
residuals are large. With;^ 2 = 2.81 (dof = 4) and thus a null- 
hypothesis probability of 2.4%, this fit can hardly be consid- 
ered statistically acceptable. However, other models, like a sin- 
gle power law (x 2 r = 7.32, dof = 6) or a bremsstrahlung model 
Of 2 = 4.77, dof = 7), result in significantly worse fits. This is 
due to the very low flux in the 40-50 keV bin with respect to 
the 20-40 keV spectrum, requiring a cutoff. 

A bremsstrahlung fit over the 20-40 keV interval yields 
a temperature of (13.6 ± 1.5) keV with^ = 1.44 (dof = 5), 
which is (within errors ) comparable to the temperature found 
bv lOwens etd] fi*999). However, over the 20-50 keV energy 
range the best-fit temperature is 9.8 ± 0.5 keV. This tem per- 
ature is in agreem ent with reports by iKubo et al.l 41998) and 
Smith et al. ( 1998), but the quality of the bremsstrahlung fit to 
our spectrum is unacceptable. 

The spectral properties of y Cas are more reminiscent of a 
typical HMXB spectrum than a WD-binary spectrum, leaving 
the question about the nature of the compact companion still 
open. 



5.3. RXJ01 46.9+61 21 

RX J0146.9+6121 (V* V831 Cas), a BeXB located 24' from 
AXP 4U 0142+61 ( see section 14.11 . was discovered 
by iMotch et aD Jl99ll) with ROSAT. Its companion star 
(LS I +61 235 ) was spectrally cl assified as a BlVe by 
iReig et alJ Jl997l) . IWhite et alJJl987l) discovered a 25 minutes 
periodicit y which was th en attributed to AXP 4U 0142+61. 
However, Hellierl i 19941) proved that this periodicity origi- 
nated from RX J0146.9+6121, using the ROSAT-PSPC imager. 
RX J0146. 9+612 also shows long-term, up to an order of mag- 
nitude, variability llHaberl et alJl99 8a). 

X -ray spectra derived from ASCA, RXTE and B eppoSAX 
data l lHaberl et alJll998blat iMereghetti etai]l2000l) could be 
modelled up to ~25 keV with power laws with or without 
high-energy cutoffs or with bremsstrahlung models. However, 
the spectra derived from the non-imaging instruments onboard 
RXTE and BeppoSAX are contaminated by 4U 0142+61. 
For the soft X-ray spectrum (<10 keV) the contribution of 
4U 0142+61 was corrected for by including fixed model pa- 
rameters, but for the hard X-ray spectrum (>10 keV) the con- 
tribution of 4U 0142+61 was still unknown. The GIS instru- 
ment onboard ASCA was able to resolve RX J0146. 9+6121 
and 4U 0142+61 in the 0.5-10 keV band. 



5.3.1. RXJ01 46.9+61 21: INTEGRAL results 

Like AXP 4U 0142+61, this source was always in the partially- 
coded FOV of IBIS-ISGRI with a total exposure of 1.13 Ms. 
With the high spatial resolution of INTEGRAL it is possible 
to resolve RX J0146.9+6121 and the AXP, and therefore there 
is no source confusion. In the 20-50 keV mosaic the source 
is detected with a detection significance of 5.3<x with a corre- 
sponding photon flux of (1.1 ± 0.2) x 10" 5 ph cm ^s^'keV -1 . 
In narrow-band (5 keV) mosaics the source is only detected 
in the 20-25 keV and 25-30 keV bands with photon fluxes 
(3.5 + 0.7) x 10~ 5 ph cirrVkeV" 1 (4.8cr) and (1.8±0.6)x 10~ 5 
ph crrr 2 s _1 keV _1 (3.3cr), respectively. In the 30^40 keV band 
RX J0146.9+6121 is still (just) visible with (0.8 + 0.3) x 10~ 5 
ph cirrVkeV" 1 {2.1(f), but in the 40-50 keV band only a 2cr 
upper limit of 4.2 x lO^ph cm 2 s ^eV 1 could be extracted. 
We fitted a power-law model to these data points to obtain an 
indication for the photon index and flux. The photon index is 
r = 3.4 ± 0.9 with a flux (1 .5 + 0.2) x 10" 1 1 erg cirT V 1 (20-40 
keV). 

Comparing the se results with the previously published re - 
sults from lHaberl etall dl998blah and lMereghetti et alJ ll200Cl) . 
our data points nicely extend the previously observed spec- 
tra. Leaving the normalization as the only free parameter, 
acceptable fits can be found for all reported fits, except 
for one, namely the single power -law fit with photon index 
T = 1.46 observed with ASCA l lHaberl et alJll998bh . does 
not give an acceptable fit to our data points. However, the 
bremsstrahlung model fitted to the same ASCA data returns 
the best slope for the INTEGRAL energy band between 20 
and 40 keV. The fluxes in the 0.5-10 keV band are also 
comparable with ~3.9 x 10 -11 erg cm^s -1 for ASCA and 
(2.9+0.6)xlO n ergcnrV 1 for INTEGRAL (extrapolated to- 
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* 
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Fig. 10. Series of 20-50 keV mosaics centered on the coordinates of IGR J00370+6122. The source was discovered in the 
revolution- 147 mosaic as a 12.3cr source. In each panel the corresponding orbital phases of each revolution are given. Both 
detections fall within the expected epoch of maximum flux. 



wards lower energies). The power-law with a high-energy cut- 
off model (r = 2.05, £„ ltnff = 17.2, E fo id = 7), fitted to RXTE 
data iMereghetti et al . (2000), gives a slightly too steep spec- 
tral slope in the INTEGRAL energy range. However, their nor- 
malization of this power law is in full agreement with the to- 
wards lower energies extrapolated INTEGRAL normalization 
of (1.4 ±0.2) x 10~ 10 ergcrrrV 1 . 

Recently, Fili ppova et alJ J2005I) reported the detection of 
RX J0146.9+6121 with independent INTEGRAL observa- 
tions. They reported a weak detection of ~3 mCrab (18-60 
keV) and a spectrum fitted with a single power law with photon 
index T = 2.9*H. The slopes of our spectra are fully consistent, 
however during their observation the source was about twice as 
bright. 

5.4. IGRJ00370+6122 

5.4.1. IGR J00370+6122: Discovery and identification 

IGR J00370+6122 is the first newly found INTEG RAL source 
in the Cassiopeia region Jden Hartog etalJE0 04b). The source 
did not show up in the integrated mosaics of revolutions 142- 
148 (1.6 Ms), but it did appear in two single revolutions. For 
this paper the data were reanalyzed with OSA 4.2, confirming 
our initial results (see Fig. 1101. The source is significantly de- 
tected in revolution 147 with a significance of 12.3cr and a flux 
of (3.1 ±0.3) x 10~ 5 ph cirrVkeV" 1 (4.6 mCrab) in the 20-50 
keV energy band. In the other revolutions the source is not de- 
tected except in revolution 142 where it appeared more weakly 



with a significance of 5.6<x and a flux of (1.4 ± 0.3) x 10 
ph crrrVkeV- 1 (2.1 mCrab, 20-50 keV). 

We redetermined the source position in the 20-50 keV mo- 
saic of revolution 147: R.A. = h 37 ra 07 s .3, Decl. = +61°21'50" 
(J2000) with an estimated 2' accuracy (90% confidence), 
which is Q'37 from th e originally distributed coordinates 
dden Hartog et al!2004bl) . 

A search through SIMBAD yields a p ossible ROSAT X - 
ray counterpart, 1RXS J003709.6+612131 JVoges et alll999l) . 
located only 0'43 from the updated coordinates. The ROSAT 
source has an optical counterpart within its 17" error radius, 
namely a ~1 0th-magnitude (B and V bands) Bllb supergiant, 
BD +60 7 3 teutledgf etalJbOOOl) . Follow-up optical obser- 
vations by iReig et all (2005) reassessed the spectral type of 
BD +60 73 to be BN0.5II-III, which belongs to the bright giant 
class. The estimated distance to the source is ~3 kpc. No radio 
emis sion is found from this source in the NVSS online archive. 
Also lCameron et ahllEcX M) did not find a radio signal with new 
VLA observations. 

We obtained a RXTE All-Sky Monitor (RXTE-ASM) light 
curve (MJD 50087-53124) and found a periodicity of (15.665 
± 0.006) days. The folded light curve in Fig. II llshows a single 
peak with a peak flux of 3.3 mCrab and a quiescent level below 
1 mCrab (after subtracting a 1 mCrab systematic bias). The 
peak flux remains above 2 mCrab for ~3 days. The epoch of 
maximum flux during the INTEGRAL and contemporaneous 
RXTE-ASM observations is MJD 53001.7 ± 0.3. A HMXB 
with an eccentric orbit where the compact object is accreting 
from wind accretion near periastron, can exhibit a light curve 



P.R. den Hartog et al.: INTEGRAL survey of the Cassiopeia region in hard X rays 



11 



0.4 




0.0 [ i i i i I i i i i I i i i i I i i i i ] 

0.0 0.5 1.0 1.5 2.0 

Phase 

Fig. 11. Folded RXTE-ASM light curve (1.1-12 keV) for 
IGR J00370+6122, including a =s 0.07 cts s -1 systematic bias, 
showing a quiescent level below 1 mCrab and a peak flux of 
3.3 mCrab (^0.23 cts s _1 ). The source is above its 75% peak 
flux only for 15% of the time, a time window of ~2.5 days per 
15.665-day cycle. Two periods are shown for clarity. 

like in Fig. ^2 Therefore this periodicity can be interpreted as 
the orbital period of a binary system. 

Using this ephemeris we found that the detections during 
revolutions 142 and 147 fall within orbital phases 0.97-1.12 
and 0.91-1.07, respectively, both at epochs of expected max- 
imum flux. This shows that IGR J00370+6122 is indeed the 
hard X-ray counterpart of the variable source which was found 
with RXTE-ASM and is most likely 1RXS J003709.6+6 12 1 3 1 . 

To derive an IBIS-ISGRI spectrum we selected only ScWs 
falling within the orbital-phase window 0.8-1.2, resulting in 
489 ks of exposure. IGR J00370+6122 is detected up to 80 
keV. The spectrum is shown in Fig.^l an d can be fitted with a 
single power-law model with a photon index of F = 1.9 ± 0.3 
(see Tabled. 

5.4.2. IGR J00370+6122: BeppoSAX archival 
observations 

BeppoSAX jBoella et alJl 1997b . operational for six years be- 
tween April 1996 and April 2002 , observed th is field with the 
Wide-Field Cameras (WFCs; ljager et all 19971) serendipitously 
on a regular basis, resulting in a total exposure of 1.37 Ms. We 
revisited these data and detected IGR J00370+6122 on four oc- 
casions in one-day averaged time bins. The spectral data were 
extracted in 31 channels between 2 and 28 keV. Channels 1, 
2 and 31 are not well calibrated and were excluded. The 2- 
28 keV fluxes of the detections are 3.0 + 0.7 mCrab (MJD 
50840. 1776-50841 .3 167), 8 .7 + 1 .4 mCrab (MJD 50980. 1576- 
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Fig. 12. Spectra (including absorption) in an E 2 F representa- 
tion of IGR J00370+6122 measured with three instruments 
at three different times: BeppoSAX Wide-Field Camera (2— 
17 keV, upper spectrum); RXTE PCA (3-20 keV, lower spec- 
trum); and INTEGRAL IBIS-ISGRI (20-80 keV). For the 
time-averaged IBIS-ISGRI spectrum only data within the or- 
bital phase bin 0.8-1.2 was used. The RXTE and BeppoSAX 
spectra were taken around phases 0.91 and 0.02, respectively. 
Overplotted are the best fits to the data which are described in 
Table [6] The order-of-magnitude flux difference between the 
RXTE and BeppoSAX observations clearly shows the variable 
nature of the peak luminosity. 



5098 1 .2459), 14.9 ± 2 . 1 mCrab (51135.81 08-5 1 1 36.4929) and 
4.6 + 0.5 mCrab (MJD 51371.9518-51373.0377). All four de- 
tections are within 0. 1 orbital period from the expected peak 
as derived from the ASM light curve. The WFCs have cov- 
ered 23 more epochs between orbital phase 0.9 and 1.1, but the 
source was not detected in any of these observations with upper 
limits less than 2 mCrab. This indicates the variable nature of 
the peak maxima. The spectrum of the brightest maximum is 
shown in Fig. El The data are fitted with an absorbed power 
law (F = 2.6 + 0.5) (see Tableg}. 

5.4.3. IGR J00370+6122: RXTE ToO observations 

In order to search for a possible X- ray pulsar we pr oposed a 
public ToO observation with RXTE ( Bradt et al.lll993l) . which 
was approved. Two short observations were performed on 
2004-05-19 and 2004-06-02 (928 s and 3728 s; observation 
ID: 90415-01). The first observation was performed near or- 
bital phase 0.14, which is during the decaying part of the 
peak. The second observation was performed near orbital phase 
0.02, which is at the maximum of the peak. In both observa- 
tions, however, the source was found to be in a low state with 
fluxes of 1.4 and 1.1 mCrab, respectively. In Fig.^Jthe back- 
ground subtracted spectrum of the longest RXTE observation 
is shown. The data are modeled with an absorbed power-law 
(F = 2.74 + 0.18, see Table|6j. Within the errors the spectra of 
both observations are similar. 
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Table 6. Spectral fit results for IGR J00370+6 1 22 -f ■ 



INTEGRAL: Power law (20-80 keV) 
T 1.9 ±0.3 

Fpoweriaw (4.0 ± 0.4) x 10-" erg crrrV 

^(dof) 1.62(3) 

BeppoSAX: Power law (2-20 keV) 
6) x 10 22 cnT 2 



N H 
I 



(13: 
2.6: 



:0.5 

** 'power law 

(30±7)x 10-" 
X 1 (dof) 1.00 (21) 



erg cm 2 s 1 



RXTE: Power law (3-20 keV) 



N H 

r 



(6.9: 

2.9i 



- 1.8) x 10 22 cnr 
0.2 



^"power law ( 1 ■ ^4 

X 2 r (dof) 0.45 (39) 



0.13) x 10"" ergcirrV 1 



5.4.4. IGR J00370+6122: Variability and absorption 

The INTEGRAL, BeppoSAX and RXTE spectra in Fig. [H] 
clearly show that the activity of the source varies widely 
from orbit to orbit. We compare the three observations in the 
5-20 keV band, to minimize effects due to possible absorption 
variations. During the (brightest) BeppoSAX observation, the 
flux was (10.7 + 1.7) x 10 -11 erg cm _2 s _1 , while during the 
RXTE observation, the flux was an order of magnitude less 
with (9.2 ± 0.6) x 10~ 12 erg cm^s -1 . The (towards lower en- 
ergies extrapolated) flux estimated from the INTEGRAL ob- 
servations falls in between the BeppoSAX and RXTE observa- 



tions with a flux of (3.4 ± 1.1) x 10" 11 erg cirT 2 ^ 1 . 

Comparing the BeppoSAX and the RXTE spectra, i.e. the 
most extreme flux measurements, show similar power laws, 
r = 2.9 ± 0.2 for the RXTE data and F = 2.6 ± 0.5 for the 
BeppoSAX data. The fitted absorption (7Y H ) to the BeppoSAX 
data of (13 + 6) x 10 22 cirT 2 is a factor of two higher than the 
fitted absorption to the RXTE data of (6.9 ± 1.8) x 10 22 cirT 2 . 
However, due to the large error on the BeppoSAX value it is 
impossible to verify whether the intrinsic absorption differs 
from one maximum to another. 



5.5. V709 Cassiopeia (V709 Cas) 

V70 9 Cas was reco gnized as an Intermediate Polar (IP) 

by lHaberl & Motchl d 19951) . IPs represent a subclass of 

Cataclysmic Variables (CVs) which consist of a weakly 

(B < 10 6 G) magnetized White Dwarf (WD) showing 

asynchronous rotation and accretion from a Roche-lobe filling 

I II I 

low-mass companion star (see Patterson 1994, for a review). 

lHaberl & Motchl i 19951) found the WD spin period of 312.8 
sec onds. The orbital per i od wa s determined to be 5.34 hours 
by iBonnet-Bidaud et all J200ll) . IPs are hard X-ray sources. 
The hard X-ray photons are believed to originate from the post- 
shock region above the magnetic poles. Their X-ray spectra can 
be described by a thermal optically-thin bremsstrahlung model 
with p lasma temperatures of tens of keV. Ide Martino et alJ 
<200 ill observed V709 Cas with BeppoSAX and RXTE. The 
time-averaged spectrum can be modeled with a plasma of 
temperature kT = 27 keV, plus a reflection component from 
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Fig. 13. IBIS-ISGRI spectrum of V709 Cas in an E 2 F represen- 
tation. Over plotted is the best fit to the data, which is a single 
power-law (see Table 0. In the bottom panel the residuals of 
the fit are shown. 

Table 7. Fit results for V709 Cassiopeia 



Power law (20-70 keV) 


r 


3.00 ±0.12 




^power law 


(4.62 ± 0.18) x 10 


-11 -2—1 

erg cm -s 


xl (dof) 


2.18(8) 




Bremsstrahlung (20-70 keV) 


kT 


19.4 ± 1.6 keV 




norm 


(2.0 ± 0.2) x 10~ 2 




(d0f) 


2.49 (8) 






(4.47 ± 0.19) x 10 


-11 -2—1 

erg cm -s 


Bremsstrahlung (20-50 keV) 


kT 


17.1 ± 1.5 keV 




norm 


(2.4 ± 0.3) x 10- 2 




X 2 r (d0f) 


1.37 (6) 




*tinabs 


(3.82^) x 10-" 


erg cm 2 s 1 



For all models the Galactic absorption was fixed at 
0.46 x 10 22 cnr 2 



the white-dwarf surface. The BeppoSAX PDS could observe 
V709 Cas up to ~70 keV. The spin modulation was observed 
up to ~25 keV with a decreasing amplitude at higher energies, 
which is typical for IPs. 

5.5.1. V709 Cas: INTEGRAL results 

With a total exposure time of 1.6 Ms, V709 Cas is detected up 
to ~70 keV (see Fig.1131. The best achieved fit of this spectrum 
is a power-law model (F = 3.00 ± 0.12, see Tabl^}. However, 
with^f 2 =2.18 (dof = 8) the fit is not really acceptable. When 
fitted with a thermal-bremsstrahlung model the fit is similarly 
unsatisfactory with x 2 r — 2.49 (dof = 8). The high x]- is due 
to the flux values in the 40-50 keV and 50-60 keV energy in- 
tervals, but there is no reason to discard these, since we did 
not find similar discrepancies in our sample of source spec- 
tra. Nevertheless, an acceptable fit can only be achieved when 
the energy bins above 50 keV are ignored. The spectrum of 
V709 Cas in the 20-50 keV energy range can then be described 
by a thermal-bremsstrahlung model (xl - 1 -37; dof = 6), with a 
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plasma temperature kT — 17.1 ± 1.5 keV. However, the plasma 
temperatures for the bremsstrahlung fits in bot h energy ra nges 
are co mparable. Compared with the results of de Martino et all 
J200ll) . the temperatures are significantly lower. However, they 
used a bremsstrahlung model including a reflection component 
from the WD. A reflection bump and a higher plasma temper- 
ature can mimic the spectral slope of a simple bremsstrahlung 
model with a lower plasma temperature. Fittin g our spectrum 
with the spectral shape modeled bv Ide Martino et all J2001 ) 
(model SAX #2; kT = keV including a reflection compo- 
nent), leaving the normalization as only free parameter, gives 
an unacceptable fit with xj- = 3.4 (dof = 9). The towards 
lower energies extrapolated unabsorbed flux (2-10 keV) of 
(3.92 + 0.15) x 10~ n erg cirrV 1 is similar (-12% lower) to 
the flux measured with BeppoSAX. Adding the plasma tem- 
perature as second free parameter, the fitted temperature is 
20.5 ± 1 .7 keV (somewhat higher than without a reflection com- 
ponent) with an unabsorbed flux (5.4±0.7)x 10~ n erg cirT 2 s _1 
{ X 2 r =2.8; dof =8). 

Fala nga et alJ J2005al) . in a parallel study of a shorter IBIS- 
ISGRI observation of V709 Cas (revolutions 262 and 263 for 
181.9 ks), derived a spectrum which they can also describe 
with a bremsstrahlung model. Their plasma temperature {kT = 
25.5^gj keV) is higher, but statistically consistent with our re- 
sults. 



5.6. IGRJ00234+6141 

IGR J00234+6141, the second new INTEGRAL source de- 
tected in the Cassiopeia field, has no proper identification yet. 
The source was detected in the mosaics for all available data, 
with a total exposure of 1.6 Ms. IGR J00234+6141 is seen as 
a weak source (0.72 + 0.12 mCrab) with a significance of 5.9cr 
in the energy band 20-50 keV and with a significance of 4.1 cr 
(1.4 + 0.3 mCrab) in the 50-100 keV band, corresponding to 
(4.8 + 0.8) x 10~ 6 ph cirrV'keV- 1 and (1.7 ± 0.4) x 10~ 6 
ph cm _2 s _1 keV~ 1 , respectively. The best source position is 
R.A. = h 23 m 24 s , Decl. = +61°41'32" (J2000) with an esti- 
mated 3' accuracy (90% confidence level), derived in the 20- 
30 keV mosaic where the source is detected with 7.4cr signifi- 
cance. 

Archival searches through SIMBAD yielded a possible 
ROSAT soft X-ray counterpart, 1RXS J002258.3+6141 1 1 
JVoges et alJ Il999t) . which is located 3.' 15 from the 
INTEGRAL centroid. It is marginally consistent with the 
INTEGRAL position. The RXTE-ASM revealed a marginal 
detection, but no modulation was found in the light curve. 
No reliable optical or IR counterpart can be found in the 
INTEGRAL error box, because of the crowded field. Online 
available radio data from NVSS showed no radio source. 

Given its position close to the Galactic plane (b = -1?0), the 
lack of a radio counterpart and the nature of the majority of 
previously reported new INTEGRAL sources near the Galactic 
plane, it is more likely to be an X-ray binary, than an extra- 
galactic source. 



Table 8. Fluxes and upper limits (UL) of weakly and non- 
detected sources. 



Source (Class): Exposure 




Energy band Flux (ph cm _2 s _1 keV~ 1 ) 


Significance 


4U 2206+54 (BeXB): 731 ks 




20-30 keV (3.6 ± 0.8) x 10~ 5 


4.7(7- 


30-^0 keV (1.2 ± 0.5) x 10" 5 


2.4cr 


1ES 0033+595 (BL-Lac): 1.6 Ms 




20-30 keV (1.3 ± 0.2) x 10^ 5 


6.3tr 


30^10 keV 2.7 x 10~ 6 


2a- UL 


4U 01 15+63 (BeXB): 1.58 Ms 




20-50 keV 2.7 x 10~ 6 


2a- UL 


50-100 keV 1.2 xlO" 6 


2a- UL 


100-150 keV 9.6 x 10~ 7 


2a- UL 


IGR J00291+5434 (accreting msec pulsar): 1.6 Ms 


20-50 keV 1.7 xlO" 6 


2a- UL 


50-100 keV 4.3 x 10" 7 


2a- UL 


100-150 keV 7.2 x 10~ 7 


2a- UL 


IGR J01363+6610 (BeXB): 1.1 Ms 


20-50 keV (5.6 ± 1.5) x lO" 6 


3.7a- (a) 


50-100 keV 1.4 x 10~ 6 


2a- UL 


<n) likely artefact 



6. Weakly and non-detected sources 

In this section information can be found on additional sources 
that have been weakly detected during this Cassiopeia-region 
survey, or have been detected with INTEGRAL at other epochs. 
The sources 4U 2206+54 and 1ES 0033+595 are only de- 
tected in one energy band with a significance greater than 5<x. 
4U 0115+63, IGR J00291+5434 and IGR J01363+6610 were 
not detected during our observations. 

6.1. 4U 2206+54 

BeXB 4U 2206+54 is a hard X-ray source (Gia cconi et alJ 
Il972i) with BP +53 2790 as optical counterpart 
(Ne gueruela & Reigl l200ll) . The orbital period is 9.6 days 
llCorbet & Peeld l200lh . Observational evidence indicate 
that the compact companion is most l ikely a neutron star 
JCorbet & Peelel l200lt iTorreion et aD l2004 iMasetti et all 
2004: lBlav et all2005h 

4U 2206+54 is located more than 11 degrees from 
Cassiopeia A. Therefore, the available exposure of 73 1 ks was 
completely obtained in the partially-coded FOV of the IBIS- 
ISGRI detector. Nevertheless, INTEGRAL has detected this 
source with a significance of 5.7<x and a flux of (1.8+0.3)x 10 -5 
ph cirrVkeV -1 (20-50 keV; see Tables|2]and|8). For a more 
detailed study of 4U 2206+54 with INTEGRAL, we refer to 
iBlav et al.M2005l) . 

6.2. 1ES 0033+595 

1ES 0033+595, an AGN of the BL- Lac type, was detected a s 
a hard X-ray source by BeppoSAX dCostamante et alJl200 lh. 
This s ource was also found with INTEGRAL bv lKuipereUdl 
d2005l) with a photon flux of ( 1 . 1 + .2) x 1 0~ 5 ph cnT 2 s" 1 ke V" 1 
(20^15 keV). We detected it in the 20-50 keV band, with 
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5.2o" significance (see Tables [5] and 03 fo r our flux values) . 
Comparing the hi gh-ene rgy fluxes of ICostamante et alJ (J2001) 
and iKuioer et all d2005b with our flux, shows the variable na- 
ture of this BL Lac. While the s ource was a factor of 2.4 
brighter during the observations of IKuioer et alJ J2005b. it was 
4.7 tim es brighter during the observations of Costam ante et alJ 
fcOO lb with a photon flux of 2.1 x 10 5 ph cirT V^eV" 1 (20-50 
keV) in December 1999. 

6.3. Non detected 4U 01 15+63, IGR J00291+5434 
and IGR J01 363+6610 

In this field three other sources were detected by INTEGRAL 
at other epochs. 4U 0115+63 was in outburst late Au gust 
and S eptember 2004 jLutovinov et alJ2004llZurita Heras et"al] 
12004b . reaching a peak flux of about 400 mCr a b. Fo r 
INTEGRAL-ToO results we refer to iFerringo et all §006). 
IGR J0029 1+543 4 was discovered by INTEGRAL early 
December 2004 dEckert et alJ l2004b . This new transient 
turned out to be the fast est accretion-powered X-ray pul- 
sar known to date (Ma rkwardt etalJ l2004blab . More de- 
tailed INTEGRAL and RXTE results are published by 
IShaw et alJ J2005b : iFaTanga et al.1 J2005bb and iGallowav et alJ 
J2005b . IGR JO 1 363 +66 10 is a short tran sient discovered with 
INTEGRAL by iGrebenev et ai] J2004 in April 2004. The 
source was seen only once for 2.3 hours. We have not detected 
any of these sources neither in the total exposure mosaics, nor 
in single revolutions, nor in single SCWs. For flux upper limits 
see Table 00 

7. Summary 

We report results from a first deep observation of the 
Cassiopeia region in the hard X-ray energy window of 20- 
200 keV, exploiting the good imaging capabilities and the 
wide field of view of the IBIS-ISGRI telescope aboard 
INTEGRAL. We report results on nine detected sources of vari- 
ous kinds: Anomalous X-ray Pulsar 4U 0142+61, several high- 
mass X-ray binaries, intermediate polar V709 Cas, BL-Lac 
IE 0033+595, and a still unidentified new INTEGRAL source 
IGR J00234+6141. In addition the upper limits are presented of 
the non-detected AXP IE 2259+586 and three more high-mass 
X-ray binaries. 

The discovery in this work of the very hard spectral tail 
between 20 and 150 keV of AXP 4U 0142+61, increases the 
number of IBIS-IS GRI AXP detections to three, the other s be- 
ing IE 1841-045 (Mo lkoy et alJl2004lKuir)er et alJl2004b and 
1RXS J170849-400910 jRevnivtsev et alJl2004X This sug- 
gests that such a hard, highly luminous spectral component is a 
common property of AXPs. Of these three AXPs, 4U 0142+61 
appears to exhibit the most extreme spectral shape (see Fig. [3}: 
below 10 keV it has the softest spectrum with a black-body 
temperature of kT = 0.470 + 0.008 keV and a very soft power- 
law component with photon index 3.40 ± 0.06. Above 10 keV 
this drastically changes into a very hard power-law compo- 
nent with photon index 0.73+0.17 (possibly even as hard as 
-0.1 1 ± 0.43 till ~70 keV if the indication for a break around 
that energy is real, see the fits in Fig. [3}- The X-ray lumi- 



nosity between 20 and 100 keV for the single power-law fit 
is 5.9 x 10 34 erg s , a factor of 440 higher than the rota- 
tional energy loss, amounting 1.3 x 10 32 erg s _1 (assuming 
M = 1.4M , R = 10 km for a homogeneous sphere and a 
distance of 3 kpc). Such a high hard X-ray luminosity has to 
find its origin in the magnetosphere of the magnetar, fed from 
the magnetic energy reservoir in the strong magnetic field. The 
underlying scenario of production processes is not yet well 
understood. To constrain different model scenarios, more de- 
tailed observational parameters are needed. One of these is the 
detailed spectral shape of the hard X-ray component. The re- 
ported COMPTEL upper limits for energies between 0.75 Me V 
and 30 MeV already indicate that the power-law spectrum of 
4U 0142+61 has to bend/break below this energy window. So 
far AXP IE 2259+586 has not been detected with INTEGRAL, 
however, the flux upper limits extracted from INTEGRAL and 
COMPTEL data do not exclude the existence of a high-energy 
tail. 

INTEGRAL has increased the number of known super- 
giant X-ray binaries s i gnificantly (see e.g. iKuulkersI [2005 ; 
Neg ueruela et al J 12006: Lut ovinov et al J 1200 61 We introduce 
a new SXB, i.e. IGR J00370+6122. We detected this source 
in two INTEGRAL revolutions, at a position consistent with 
that of the ROSAT source 1RXS J003709.6+612131, which 
has a BN0.5II-III bright-giant counterpart. Analysis of archival 
RXTE-ASM data revealed the 15.665 ± 0.006 day orbital pe- 
riod. Follow-up studies using a ToO RXTE observation and 
BeppoS AX-WFC archival data rendered a consistent picture of 
a new, by an order of magnitude variable, wind-accreting SXB. 
The source exhibits below 20 keV a power-law spectrum with 
index ~2.7, and a hard power-law spectrum with index 1.9 + 0.3 
in the INTEGRAL window 20-80 keV. 

We report the detection of a second SXB, 2S 0114+650. 
For the first time the spectrum has been measured up to ~150 
keV. The good statistics in the time-averaged spectrum above 
70 keV, showed that a single bremsstrahlung model cannot fit 
the data satisfactorily. In addition, a power-law component with 
index 2.0 ± 0.7 is required, or the spectrum can be explained 
with a single power law with index 3.01 ± 0.06. 

Our deep survey of the Cassiopeia region rendered also de- 
tections of three Be X-ray binaries: y Cas, RX J0146.9+6121 
and 4U 2206+54. For the first two, the INTEGRAL spectra 
were measured to higher energies than published earlier, and 
are discussed in comparison to these earlier reports. 

The intermediate polar V709 Cas was detected up to ~70 
keV. The best model fit was obtained for a power-law model 
yielding a photon index of 3.00 + 0.12, different from the 
bremsstrahlung-model fits reported in literature. 

Finally, two more hard-X-ray sources are found: 1) 
IGR J00234+6141, discovered as a weak (~1 mCrab) source 
up to 100 keV. No counterpart could be found yet. Since it is 
located close to the Galactic plane and there is no radio coun- 
terpart, it is likely an X-ray binary, rather than an AGN. 2) 
1ES 0033+595, an AGN of the BL Lac type, also detected in 
this work at the 1-mCrab level. 
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